


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1972-03 


The measurement and correlation of sound 
velocity and temperature fluctuations near 
the sea surface 


Duchock, Charles Jack 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/16249 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
A8 DUDLEY research materials and institutional publications created by the NPS community. 
dil əra Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


il KNOX appointed — and published — scholarly author. 

ğa. LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


THE MEASUREMENT AND CORRELATION OF 
SOUND VELOCITY AND TEMPERATURE 
FLUCTUATIONS NEAR THE SEA SURFACE 


Charles Jack Duchock 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


THE MEASUREMENT AND CORRELATION OF SOUND 


VELOCITY AND TEMPERATURE FLUCTUATIONS 


NEAR THE SEA SURFACE 


by 


Charles Jack Duchock, Jr. 





March 1972 


Approved gor public release; distribution unlimited. 


a 





The Measurement and Correlation of Sound Velocity 


and Temperature Fluctuations Near the Sea Surface 


by 


Charles Jack Duchock, Jr. 
Lieutenant Commander, United States Navy 
B.S., Auburn University, 1964 


Submitted in partial fulfillment of 
requirements for the degree of 


MASTER OF SCIENCE IN OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1972 





ABS TRACT 


The temporal variation of salinity, temperature, water 
turbulence components, surface wave height, acoustic ampli- 
tude and acoustic velocity were studied statistically by 
computing the auto-correlations and power spectral densities. 
Correlation times were high for temperature fluctuations 
with a decay time (el) being of the order of forty seconds 
or greater for temperature and of the order of two seconds 
for the other parameters. All parameters were shown to 
exhibit maximum energy at the periods of the predominant 


surface wave energy. 
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I. INTRODUCTION 


This work describes part of a joint experiment by 
members of the Physics and Oceanography Departments of the 
Naval Postgraduate School. The Oceanographers basically 
described the environmental parameters and fluctuations in 
them while the Physicists simultaneously recorded the affect 


of the environment on sound projected into it. 


A. OBJECTIVES 

The major objectives of this work were to describe the 
oceanic parameters and fluctuations in them as completely as 
possible while simultaneously recording fluctuations of a 


sound signal projected into the water. 


1. Knowledge of the Environment 


Whereas we seek to understand the complex inter- 
relationships of the air and ocean, of particular importance 
is the affect of the oceanic environment on tnose devices 
and energy sources put into it by man for his various 
purposes. The particular case to be investigated here is 
the affect of the ocean on sound and the interrelationships 
Of sound and certain environmental parameters near the 
surface. 

a. Oceanographic Environmental Factors Affecting 


Sound 


Traditionally, salinity, temperature, and 


pressure are the parameters considered as affecting the 
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velocity of sound in sea water. A more comprehensive 
knowledge of the environment into which we project sound 
could show that other factors have a significant role in 
determining the velocity of sound. 

Near the ocean surface the salinity is changed rapidly 
as evaporation, condensation, and precipitation occur. The 
affect of the wind on the surface coupled with changes in 
density cause mixing of the water to occur at varying rates 
depending on the differences in density in the water column 
and the velocity of the wind as well as its duration and 
fetch. This mixing is, however, never so thorough as to 
cause the water to be homogeneous in the near surface layer. 

The temperature of the surface layer can be changed by 
insolation. The amount of heat added to the water by insola- 
tion is dependent upon the time of day, cloud cover and the 
elevation of the sun. The temperature of the water at the 
air-water interface is constantly changed by heat transferred 
between the air and water masses. The direct heat transfer 
is a function of the difference in temperature that exists 
between the air and the water. An air mass warmer than the 
water will transfer heat to the water at the interface and 
exert a stabilizing influence on the water column. An 
overlying cold air mass has the opposite affect and tends to 
Cause mixing of the water by convective overturn. Wave 
action also aids the mixing process but again the water is 
not usually mixed to the extent that it becomes homogeneous 


but is found to be in a constant state of flux due to these 


processes. 
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It has long been known that in the ocean medium sound 
pulses propagating over relatively short paths are phase and 
amplitude modulated due to scattering from inhomogeneities 
such as the temperature microstructure in the ocean. These 
thermal inhomogeneities produce random variations in the 
refractive index of the medium. Thus there is a need for 
determination of the spatial and temporal variation in the 
velocity of sound with temperature and in fact with all 
parameters that affects its propagation. 

b. Interrelationship of Spectra 

A study of interrelations between spectra of 
Ocean surface heights, water particle velocities, tempera- 
pure and salinity fluctuations, acoustic amplitude fluctua- 
tions, and acoustic phase modulation could give some insight 
into the affects of inhomogeneities in the environment on 
sound.  Time-series measurements are needed to develop the 
Statistical relationships between wave spectra, underwater 
turbulence, temperature microstructure, bubble populations, 


and scattering of acoustic energy from the sea surface. 


This need was stated in 1960 when LaFond noted the necessity 
of simultaneously measuring during acoustic tests a number 
Of oceanographic parameters which were believed to influence 


underwater sound transmission. 


2. Scope of the Research as it Applies to Near Surface 
Problems 


a. Data Needed 
Data must be acquired to fully describe the 


environment into which sound is being projected. These will 
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allow comparison of results in a systematic fashion which 
could assist in determining specific associations heretofore 
unnoticed or undocumented. 

b. Purpose of Experiment 

The purpose of the experiment was to completely 
Geseribe the small-scale physical properties in the upper 
ocean and to determine their temporal and spatial inter- 
relationships. The experiment was designed to investigate 
the complete near-surface physical environment and its 
interaction with sound waves. 

In 1969 Campanella performed an experiment whose 
objective was to experimentally observe the relationship 
between the time autocorrelation function of the amplitude 
of sonic pulses propagated in a random medium and the time 
autocorrelation function of the temperature microstructure, 
and hence of the random refractive index of the random 
medium. His EE (performed in electrically heated 
tap water in a small laboratory tank) showed that the time 
autocorrelation of the sonic-pulse amplitudes was equivalent 
to that measured for the temperature microstructure. 

in realizing the object of determining the 
relation and affect of the oceanographic parameters on the 
acoustic parameters there exists a secondary objective to 
investigate the feasibility of using acoustic devices as 
sensors of fluctuating oceanographic parameters. 

Many spatial and time scales affect sound 


Propagation. However, these studies are oriented towards 
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achieving an understanding of small scale interactions in 
the upper ocean. Small means spatial scales of the order of 
meters to the smallest size able to be resolved by the 
instrumentation either by size or by the time constant of 
it. The largest scales will be imposed by the length of the 


record. 


7 Parameters to be Studied 
A determination of the specific parameters to be 
studied was based on instrumentation available and the 
limitations of the recording equipment. There are several 
obvious omissions such as bubble populations and distribution 
and particulate matter but the measurement of these para- 
meters will have to remain unrecorded until more appropriate 
instruments are available for open ocean studies. 
a. What Must Be Measured 
The basic parameters which must be measured to 


define the environment are the following: 


(1) Wave height (h) 


(2) Turbulence Components (u,v,w) 

(3) Temperature Microstructure (T) 
(4) Sound Velocity (c) 

(5) Salinity (S) 

(6) Acoustic Amplitude Modulation (A) 


(7) Acoustic Phase Modulation (40) 


b. Where Measurements Were Taken 


The actual experiment was conducted at the Naval 


Undersea Research and Development Center (NURDC) Oceanographic 
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platform, one mile off Mission Beach, San Diego, California. 
The tower has carts mounted on rails on three sides which 
can traverse over the vertical column from the surface to 
the ocean floor. Instruments were mounted on the most 
seaward cart and measurements were made at various depths 
in the water. 
c. When Measurements Were Made 

The measurements were taken on the twenty-first 

and twenty-second of October 1971 with as many parameters 


being recorded according to equipment limitations. 


B. REVIEW OF THE PROBLEM 

A complete review of the problem would necessitate 
coverage of the study of the ocean and sound in the ocean 
from earliest recorded time to the present day. The intent 
here is not to be historical but to show the basic content 
of past efforts and relate where knowledge of the subject 
stands today. 


l. General Review 


The study of sound in the ocean was initiated by 
Leonardo Da Vinci in 1490. He described a method of ship 
location using a hollow tube with one end placed beneath the 
surface of the ocean and with the human ear placed at the 
other end to allow another ship to be heard at great distance. 

a. Historical Background 

In the year 1827 Colladon and Sturm actually 


measured the velocity of sound in Lake Geneva. Many other 
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scientists worked on the study of sound in the ocean and the 
factors affecting it through the years. It became apparent 
that the major factor controlling the velocity of sound in 
water was the temperature of the medium. This fact led 
LaFond to state in 1962 that the principal oceanographic 
study to be undertaken in relation to sound velocity was the 
general temperature structural variation of the ocean with 
time and space. His pursuit of Knowledge along this line 
led to a great amount of data on internal waves in the ocean. 

LaFond's basic reason for studying internal waves 
was their affect on underwater acoustic transmission. Sound 
passing through these waves is refracted, depending on the 
Base of incidence with the thermocline containing the 
internal waves (and other irregularities) and the vertical 
temperature gradient. 

It has been determined that in the upper few 
hundred meters of the open ocean the density structure is 
almost entirely dependent on temperature variations. 
Measurements over the past few years have indicated that one 
can no longer regard temperature as a continuous variable 
with depth. Instead, extremely narrow layers (lenses) of 
water separated by strong temperature gradients appear to be 
interleaved. 

b. Most Recent Measurements and Findings 

An investigation of the most recent measurements 

and findings on the velocity of sound in the ocean and its 


relation to oceanic parameters reveals that a great deal of 
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time and effort has been spent recording sound velocity in 
various known water types and by curve fitting procedures 
developing equations to determine c knowing T, S, and P. 
If one assumes that the speed of sound changes smoothly with 
temperature and that the sound-speed derivative decreases 
monotonically as temperature increases, then aside from 
questions of absolute accuracy, the Greenspan and Tschiegg 
distilled-water equation is considered internally the most 
self consistent (Lovett, 1969). The standard deviation is 
2.6 centimeters per second over the widest temperature range 
(specifically zero to one hundred degrees centigrade). 

The exact water type being dealt with in the open 
Ocean is rarely known however, and the more that is under- 
stood about the affects of oceanic variables the better will 
be the results of attempts to predict exactly what sound 
will do when projected into the ocean. The fact that the 
ocean is far from the nearly homogeneous well mixed fluid 
that it was thought to be has been amply demonstrated by 
measurements over the past few years by Stommel and Federov 
(1967), Woods (1968) and Denner (1967). 

Since there exist many formulas for computing the 
velocity of sound in use today a look at several will be most 
helpful in a later discussion of the work carried out at the 
NURDC Oceanographic tower. Wilson's formula (Wilson, 1960) 
which was arrived at using only laboratory data from actual 


ocean S.T.P combinations is: 
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c = 1449.30+1.5848 x107l p+1.572x10-5 P2 -3.46x10712 p4 


ÜNSNIS.396 x 10-2 T2 4 2.604 x 107% q3 


+ 4.587 x LO 

+ 1.19 x 100 (s - 35) + 9.6 x 1072 (s – 35)2 

* 1.354 x 107? q?P - 7.19 х 1077 ТР? 

SDN 10-29(s"- 35) m 
Mais is accurate to six significant digits with dimensions of 
İ - meters @er second, T in degrees centigrade, S “in parts 
per thousand, and P in kilograms per square centimeter. 

Using Kattegat water with a salinity of 27.405 

parts per thousand over the temperature range 4.237 degrees 
centigrade to 25.164 degrees centigrade the least squares fit 


gave c in meters per second as: 


c 7 1438.815-* 4.69236 T4g -5.4843 x10 ^ 17,5 4 2.77 x 10 ^. T3, 
(uS 

-2 2 -4 3 
¢ =1438.813+ 4.69477 Teg -5.4926 x 10°“ T^c, 4 2.78 x10 " T". 


with an estimated accuracy of + 0.05 meters per second. 


The International Practical Temperature Scale of 
1968 was adopted by the Comite International des Poids et 
Mesures at the October 1968 meeting. The new scale replaced 
IPTS-48 which had been published by the National Bureau of 
Standards (Metrologia, 1969) as the official measure.  NBS 
Calibrations for platinum resistance thermometers. liquid-in- 
glass thermometers, and thermocouples have generally been 


based on the new scale since 1969. 
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Del Grosso's determination (Del Grosso, 1970) 
B zı nq salinities from 3l to 39 parts per thousand and 
temperatures from zero to thirty degrees centigrade is: 


с = 1449.0634 + 4.57462 T4g - 5.27147 x 107? T^,, + 


2.46419 x 1074 T3,g + (S - 35) (1.34455 - 1.32888 x 1072 ə 


* 1.0444 x 1074 T2,g). 

All of the above determinations for c emphasize 
the dependence of the velocity of sound on the temperature. 

The work done by Liebermann in 1951 showed that 
temperature inhomogeneities possessed surprisingly small 
linear dimensions in spite of diffusion processes which tend 
to obliterate them. 

Shonting found that the more closely one examines 
the temperature field in the upper layers of the ocean the 
more complex it appears. Small scale temperature fluctuations 
of the order of 0.05 to one degree centigrade occur spatially 
distributed in the surface layers. These thermal anomalies 
apparently occur as highly turbulent or eddy-like features a 
few centimeters or meters in diameter or even as thin strata- 
fied layers a few meters thick extending horizontally many 
kilometers. He used the term thermal microstructure to define 


the observations he had made. 


2. Specific Review 
Seeking information about specific open ocean studies 
of the affects of temperature on the velocity of sound we 


must return to the work done by LaFond. 
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a. San Diego Experiment Background 

LaFond realized that the temperature structure 
influenced the transmission of sound in th sea and this led 
to his studies of the vertical fluctuations of isotherms 
with respect to both time and distance. He knew that a 
knowledge of the magnitude of these fluctuations would help 
to understand if not solve the problems of sound transmission 
Variations in the ocean environment. After associating the 
vertical fluctuations of the temperature structure with 
internal waves, he found that the most significant thermal 
fluctuations, from the standpoint of equipment response and 
Variability in sound transmission, are the fluctuations 
caused by short-period internal waves. These may be grouped 
in the two minute to two hour wave period range and consist 
of rapid temperature changes with time that are caused 
largely by the vertical motion of the thermocline. 

LaFond found that the correlation of temperature with 
reference to time decreased at all depths and that in 
Semeral terms, the correlation is poor in the isothermal 
layer and significant for levels in the thermocline and 
below. He related the low correlation in the near surface 
to fluctuating atmospheric conditions and stated that this 
results in patchiness of the surface temperatures even 
though the fluctuations of the atmospheric conditions may be 


small. 
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b. What Has Been Measured 

For many analytical purposes most of the ocean 
for most of the year can be taken to consist of three strata. 
The first stratum is the surface layer consisting of well 
mixed, nearly isothermal water. The second is the lower more 
50 vater in which the temperature and density vary slowly. 
Separating these two parts is the seasonal thermocline region 
in which the vertical temperature gradient is relatively 
иә. The horizontal salinity gradient is slight and iso- 
thermal surfaces nearly coincide with isopycnal surfaces 
here and an obvious method of observation of internal waves 
Ero measure the vertical temperature or salinity structure 
in the vicinity of the thermocline.  LaFond chose to measure 
the temperature structure. 

Pederson reports simultaneous measurements of 
sound velocity, temperature and pressure along with depth 
using the self-propelled underwater research vehicle. A 
definite change in sound velocity was sensed by the sound 
velocimeter and correlated to a change in temperature of 
0.05 degrees centigrade. Apparently no attempt was made to 
perform statistical analysis of the signals recorded using 
SPURV. 

c. What Is Known About This Site 

LaFond (1962) defined internal waves as undula- 
ting swells occurring between subsurface water layers of 
varying density, even though the density change be slight. 


It is significant to note that vertical oscillations in the 
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Fu ral and density structure of the sea are apparently 
present in all oceans and at all depths. In his measure- 
ments of the vertical oscillations of the thermal structure 
LaFond found that the speed of internal waves was from 0.11 
to 0.6 knots with an average value of 0.3 knots. Fifty per 
cent of the internal waves were found to be higher than 5.6 
feet and fifty per cent of all waves longer than two minutes 
were found to have periods greater than 7.3 minutes. These 
measurements were made from the same site from which the 
data for this experiment was taken and gives some indication 
of the degree of activity at this location. However these 
measurements were made in the summer when a sharp thermo- 


cline existed. 
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II. EXPERIMENT 


The background material previously discussed provides 


the basis for a discussion of the actual experiment. 


A. SITE 
The site chosen was the Naval Undersea Research and 
Development Center Oceanographic Tower in San Diego, 
California (Figure 1). 
1. History 
The tower was installed in the summer of 1959 in 
water 18 meters in mean depth 1.6 kilometers (approximately 
One mile) off Mission Beach in San Diego, California. 
2. Description 
The tower is constructed of four tubular steel legs 
Ku (o the ocean floor vith steel pins sixty-three feet 
long driven into the sandy bottom. A pipe framevork supports 
a cement deck and an instrument house tventy-three feet above 
the vaterline. There are tvo catvalks belov deck level vhich 
are used for handling gear. Several types of instrument 
handling gear are available including rail mounted carriages 
fixed to three sides of the tower. Figure two illustrates 


the tower as originally configured for NURDC studies. 


3. Reason For Selecting 


The tower was selected for its stability, accessi- 


bility, low self-noise level, exposure to the open ocean and 
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its economy. One of the major advantages of the tower is 
that it allowed laboratory-like controls in a natural 

environment. The limitations of the tower are of course 
inherent in the fact that it is an immobile platform and 


therefore restricts measurements to relatively shallow water. 


B. SENSORS 

The sensors used in the experiment were chosen on the 
basis of availability and reliability. They were mounted so 
that the mean current flowed into the sensor array and then 


past the mounting structure. 


1. Sound Velocity Sensor 


Sound velocity was measured with a Ramsay MK-1 Deep 
Sea Probe utilizing only the sound velocimeter portion of 
the probe. 

a. Description 

Velocimeters, precision imstruments for directly 

measuring the speed of sound in water, evolved from the 
first successful sing-around prototype developed in 1957 by 
Greenspan and Tschiegg of the National Bureau of Standards. 
The Ramsay velocimeter is of the sing-around type consisting 
Of a transducer and transmitter which transmits a pulse of 
four MHz into a twenty-five centimeter sound path. The 
pulse is reflected twice to reduce errors due to water 
motion. After reflection the pulse is picked up on a receiv- 
ing transducer that is the input for a high-gain pulse- 


Shaping amplifier. The amplifier retriggers the blocking 
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oscillator, and a repetition frequency results which is 
higher than the free-running rate. Thus the water path acts 
as the delay line where the variation in sound velocity 
through the water changes the delay and hence the "sing- 
around" frequency. The frequency also depends on the path 
length and the circuit delays. The configuration of the 
path length makes it impossible to measure the length to any 
desired degree of accuracy. Also, because of selective 
attenuation, the received pulse rises slowly in comparison 
with the sent pulse and an unknown time delay is introduced 
during which the received pulse is below the noise level. 

Because of this the instrument must be 
calibrated in a liquid for which the velocity of sound is 
known accurately, and the liquid must be similar to that in 
which the instrument will be used. Thus if the instrument 
is to be used in sea water, it may be calibrated in 
distilled water in which the sound velocity is known as a 
Z netion of temperature. 

The pulse repetition frequency (PRF) is betveen 
5,600 and 6,400 Hz depending on the velocity of sound. 
Sound velocity in meters per second can be found by dividing 
50 EF by four. The PREF controls the output of the sound 
K OCity oscillator. The output of the sound velocity 
oscillator will be one half the PRF and twice the velocity 
Of sound in meters per second in the Ramsay probe used in 


this work. 
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The electronic package of the probe is contained 
within a stainless-steel pressure-proof housing that can 
withstand the pressures found at maximum ocean depths. 

The characteristics of the unit are: Range-1,400 
to 1,600 meters per second; Output Frequency - 2,800 to 
S7200 Ez; calibration accuracy —- 0.01 meters per second; 
response time - 160 microseconds. 

The signal from the probe was transmitted via 
the sea cable to the deck unit. The Ramsay Deep-Sea Probe 
Control Unit, Model 1B, amplifies the data signal, demulti- 
Plees the f-m telemetry tone signals and converts these 
ne Signals to variable D.C. voltages.  Frequency-to-DC 
converters provide zero to ten volts D.C. over the frequency 
spans of the channel bandwidths (2,800 Hz to 3,200 Hz for 
the velocimeter). A block diagram of the velocimeter system 
is shown in Figure 3. 

The signal was then filtered with a Khronhite 
Ultra-low Frequency Band Pass filter model 330-A set on high 
input amplitude with a low cut-off frequency of 0.02 Hz and 
205niuqh cut-off frequency of 2,000 Hz. The signal was then 
amplified vith a Hevlitt-Packard amplifier model HP 2470 


before it was recorded. 


b. Employment 
The velocimeter probe was mounted in the vertical 
center of the instrument array behind the transducer used 
for amplitude and phase modulation studies. This placement 


allowed measurements of sound velocity in close proximity to 


the other instruments without interfering in their operation. 
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2. Temperature 


Temperature was measured using a specially construc- 
ted Wheatstone bridge circuit utilizing a temperature sensi- 
Eve thermistor as one leg of the circuit. A total of three 
such circuits were used. 

a, Description 

Thermistor beads are small glass-encapsulated 
high-thermal-resistant materials with connecting leads. The 
resistance of the beads is greatly influenced by temperature 
(hence thermal resistor). The temperature changes the bead 
5 --nce which"echanges thé output voltage of the bridge. 

The thermistors were used as one leg of a 
Mame stone bridge circuit. These circuits were so construc- 
ted that a null coulá be set for any resistance value of the 
thermistor and therefore for any desired temperature. The 
esreuit üsed and its relation to the amplifier and recorder 
is illustrated in Figure 4. 

The temperature sensor used was a type K496 iso- 
curve thermistor manufactured by Fenwall Electronics, Inc., 
Farmington, Massachusetts. It is a double-bead type which 
incorporates a special aging process, giving stable and 
repeatable accuracies not usually obtained with ordinary 
paeuxmStors. These thermistors were taken from expendable 
bathythermographs. Statistically, if there is a drift in 
the characteristics of one bead, the other bead may not 
Change or may shift in the opposite direction. This tends 


to minimize any overall calibration change. 
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The sensor is epoxy-potted within a stainless 
steel tube three millimeters in diameter. The beads are 
coated with an extremely thin layer of glass which provides 
good electrical insulation, stops microscopic leakage of gas 
caused by slight chemical strain release in the semi- 
conductor caused by hydrostatic pressure, and most impor- 
tantly permits a thermal response time in sea water of less 
than 100 milliseconds. 

These thermistors are capable of operating over 
the range from -5 to 30 degrees centigrade making them suit- 
able for ocean temperature measurements. The thermistor"s 
relatively large resistance change per degree of temperature 
Beoviıdes a high resolution factor. It has a resistance of 
about 15,000 ohms at 25.0 degrees centigrade and exhibits a 
nominal resistance change of about 570 ohms for a one-degree 
ecntigrade change at 25 degrees centigrade as compared to 
less than two ohms for a typical 200-ohm platinum resistance 
sensor. 

The high resistance of a thermistor coupled with 
its high sensitivity makes it appropriate for use at the end 
Of a conducting cable. Changes in the resistance of the 
cable due to ambient temperature variations can be neglected 
for short cable lengths up to about 1,000 meters, since the 
resistance of the cable is small compared to the resistance 
mene thermistor. By maintaining a constant current the 
voltage across the thermistor is a direct function of the 


temperature therefore requiring no reference or cold 


29 





junctions and high amplifier gain is not required which 
reduces errors due to high gain instability. 
b. Construction 
Construction of the associated bridges vas 
accomplished with the assistance of Dana Maberry, a techni- 
cian from the Oceanography Department. 


(1) Choosing Components.  Thermistors were 


chosen müs Of their high sensitivity, fast time-response, 
minimal effects from lead resistance and lov cost. Their 
lov cost alloved several temperature sensors to be used 
simultaneously. 

(2) Assembiy. Assembly of the thermistor 
devices vas accomplished using a separate bridge box for 
Sach of three thermistors with three additional thermistors 
prepared in advance in case of breakage of the delicate 
beads. 

(3) Calibration. Calibration was done in a 
stirred Dewar flask of water using a quartz thermometer 
standard. The temperature was varied from ten degrees 
centigrade to eighteen degrees centigrade in approximately 
one-half degree increments. The circuits were nulled using 
a single decade box and the nulling resistance was plotted 
versus temperature. As the temperature increased the 
resistance of the circuit decreased. The slope of the line 
for this range was 910 ohms per degree centigrade. Figures 
five and six give typical calibration curves obtained for the 


Buerxrmistor circuits. 
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The calibration procedure was repeated with 
Mae circuit nulled at ten degrees centigrade and the voltage 
ouput recorded. The plot showed that as the temperature 
increased the voltage increased and the bridge proved to be 
very linear showing an increase of 0.0553 volts per degree 
centigrade. By use of the single decade box the circuits 
could be nulled at any desired temperature allowing fluctua- 
tions about this "baseline" temperature to be monitored. 
c. Employment 

The arrangement of instruments was chosen to 
provide for optimum coverage of the volume of water contain- 
Eug all instruments without causing mutual interference. 
The sensors were mounted as shown in Figures seven, eight 


and nine. 


3. Remaining Sensors 
The sensors described above are those which were the 
и -epomsibility of the author. The following is a brief 
description of the other sensors used in the experiment. 
a. Baylor Wave Gauge and Pressure Wave Sensor 
The surface waves (h) were measured using a 
wace penetrating resistance wave gauge manufactured by 
the Baylor Company and also with a pressure wave sensor made 
by Interstate Electronics. 
b. ®Warbulence Meter 
The turbulence measurements (u, v and w) were 
made using an electromagnetic flow meter manufactured by 


Engineering Physics Company. This instrument is capable of 
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measuring two orthogonal components simultaneously. The 
resolution of this instrument is limited by the wave number 
scale because of the volume over which the flow velocity is 
integrated. This limits this particular model to measure- 
ments of turbulence with wave lengths greater than about 
twenty centimeters. Using Taylor's hypothesis this limits 
the frequency resolution to about one Hertz. 
c. Salinometer 
A Bissett-Berman inductance salinometer was used 
to obtain the temporal variation of Salinity (S). 
d. Sound Device 
The amplitude (A) and phase modulation (@) of a 
continuous sound wave of 20-80 KHz was measured using a 
source and hydrophone separated by approximately two meters 


(reference Theses of LCDR Smith and LCDR Routmann). 


C. DATA COLLECTION 
1. Recorder Used 
All data were recorded on a fourteen channel Sangamo 
model 3500 magnetic tape recorder using FM electronics. 
Data recording was at 1-7/8 inches per second. The length 
of each run was approximately twenty minutes. 
2. Dates of Collection 
The data recording was done on 21-22 October, 1971. 
Table I is a tabulation of data collected with a check 


indicating data and a zero indicating non-collection. 
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3. Auxiliary Data Assembled 
As additional data the standard oceanographic 
observations as delineated in H.O. publication 607 were made. 
An attempt was made to observe sea life near the instrument 
array but during the hours of darkness and when the array 
was deeper than a few meters such attempts were generally 


inconclusive. 
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III. ANALYSIS OF DATA 


The analog data recorded on magnetic tape consisted of 
runs of approximately twenty minutes duration. The data 
feet Cep required transfer of this information to strip 
charts which were then used for digitization and statistical 
analysis. The strip charts themselves were useful in visual 


correlation work with the records. 


A.  DIGITIZATION 
The digitization procedures were accomplished in three 
Steps as outlined here. 


E Strip Chart Production 


The tape recorded data was collected at a tape speed 
of 1-7/8 inches per second. The tapes were replayed at 7-1/2 
inches per second (or 4x speed up) and the information 
recorded using a Brush Mark 200 eight channel strip-chart 
recorder. The increased playback Speed reduced the twenty 
minute records to five minutes in length. A speed of five 
millimeters per second was used on the Brush recorder produc- 
ing a record 1.5 meter in length. These records were used 
for visual correlation studies and were also the basis for 
Exosuction of a digitized record for statistical analysis. 


2. Digitized Magnetic Tape Production 


mme strip chart records were taken to the Fleet 


Numerical Weather Central facility at Point Pinos, California 
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where a digitized magnetic tape was produced utilizing the 
Calma Company model 480 digitizer. The digitizer recorded 
each one-hundredth of an inch variation of tthe trace from 
the starting point of each run which was marked by a verti- 


иллә imposed on the chart at the start of each recording. 


5064 CardiProğuction 
The nine track magnetic tape obtained from the Calma 
5”  “izer vas converted to data cards usinğ the CDC 5000 
Computer at Fleet Numerical Weather Central. These cards 
were then used as the input data for the statistical analysis 


done on the Naval Postgraduate School's IBM 360-67 Computer. 


B. COMPUTER PROGRAMS 

The programs: utilized to analyze the digitized data were 
from the W. R. Church Computer Center library and the files 
of Fleet Numerical Weather Center's computer center. 


l.  MISR Program 


Hie subroutine MISR from the We KR. “hurch Computer 
Center library was used to compute the mean, standard devia- 
tion, skewness and kurtosis for each of the data sets. This 
was done to serve as a comparison with the values obtained 
from the second program used in the analysis. 

2. —FNWC Spectral Analysis Program 

This program was the basis for the analysis of data 
for this experiment. The time increment in this program was 
0.2032 seconds per sample which was arrived at by considering 


the speeds of the recorders and the incremental distance 
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27 rent in the digitization process. This set the Nyquist 
frequency at 2.46 HZ for all the data analyzed. The two 
products of interest from this program were the normalized 
autocorrelation and power spectrum graphs from the CALCOMP 
plotter. 
a. Normalized Autocorrelation Graphs 

The Normalized Autocorrelations for each data 
set were plotted with an abcissa scaled in seconds to a 
Maximum time lag of forty seconds. These plots allowed a 
quick evaluation of the persistence exhibited by each para- 
meter. Generally the temperature data showed a high degree 
of coherence while all other parameters showed a much lower 
Coherence. The graphs are shown as Figures 10 to 28. 

b. Power Spectral Density Graphs 

The power spectra were plotted for each data set 
Ee show the contributions of oscillations with various 
frequencies to the variance of each time series recorded in 
the experiment. From the graphs shown as Figures 29 to 47 
it is noted that though the maximum energy is contributed by 
periods longer than one-half Hertz the total energy contri- 


Dated by periods less than this is by no means negligible. 
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IV. CONCLUSIONS 


A. STRIP-CHART COMPARISONS 

A visual investigation of the records obtained from the 
Brush recorder showed an apparently high degree of correla- 
tion between the salinity, wave profile, turbulence profiles, 
sound velocity and amplitude profiles. The periodicity of 
[ol these parameters was about that of the surface waves 
which were recorded at the same times. The signals from the 
thermistors did not appear to correlate. A possible reason 
for this observed behavior of the thermistor signals is that 
the time constant of the thermistor beads used as sensors was 
much smaller than that of the other instruments. 

in Run Three 

A detailed investigation of the records for this run 

revealed that the salinity trace and wave height were almost 
identical records. These records showed no long term fluc- 
tuations with the signals showing short period fluctuations 
about an average value. This base value for salinity was 
about 33.5 parts per thousand from the Bissett-Berman 
salinometer and 33.512 parts per thousand from analysis of 
a water Sample taken at the time of the run. This extremely 
strong correlation between salinity and wave profile at 
Shallow depths held throughout the experiment. There exists 
a strong possibility that the instrument is somewhat velocity 


dependent. The vertical orientation of the axis of the 
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cores of the inductance coils (used as the sensors) could 

have caused the water contained in the core to move only with 
the vertical component of the turbulent velocity caused by 
mavye action. This velocity may have been recorded rather 

than actual salinity fluctuations. This suspicion is substan- 
tiated by the similarities between the vertical component of 
ENS ty measured with the turbulence meter and the salinity 
and wave records as well as by the fact that the correlation 
decreased as the depth at which the recordings were taken 
increased. 

The temperature records showed a long term positive 
mecrease throughout this run and the fluctuations appeared to 
be about the same level throughout the run. The phase 
25  tude signal fluctuated vhen the temperature increased. 
This relation held for all the runs analyzed. There is some 
indication that the phase amplitude fluctuations decreased 
if the temperature remained constant in D.C. level, even 


Enough short term fluctuations continued. 


The vertical component of turbulence and the phase 
amplitude appeared strongly correlated in frequency, phase 
and amplitude at times but this correlation did not hold 
throughout the run. 

The long term changes in the temperature records 
indicates that there was a net movement of water past the 
sensors. The sound amplitude fluctuations which increased 


when there was a temperature change support this idea. 
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2. Run Four 

Tne salinity and wave records for this run show very 
little correlation. Since this particular run was made near 
the bottom and wave action here was minimal, this observa- 
tion supports the conclusion that the salinometer is velocity 
dependent. 

The apparent dependence of phase amplitude fluctua- 
tions on temperature changes is again evident in this data. 
ethe temperature showed long term changes the fluctuations 
wne phase amplitude increased. There also appeared to be 
correlation between the temperature and salinity signals for 
MTS run. The signal from the sound velocimeter also 
exhibited fluctuations at the same times as those in the 
memoerature record. 

57 ın Five 

Hn Comrelation Detween the salinity ana surface 
Mave Signals is high for this run. These data were made at 
a shallow depth as was run three, therefore supporting the 
idea that the salinometer is velocity dependent. 

mic reratcion between the fluctuations in phase 
amplitude and an increase in temperature are very evident 
once again in this record. 

fhe conclusions, or ideas, based on the visual 
observations of chart records is of course somewhat subjec- 
tive and must be considered tentative at best. However they 
point out features that are important enough to warrant 


detailed examination in future experiments. 
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B. AUTO-CORRELATION GRAPH COMPARISONS 

The most opvious fact obtained from a study of the 
auto-correlation graphs is the long correlation time exhibited 
2” ne temperature signals. Table IV lists the times for the 
signals to decay to a value of el. Тһе temperature signals 
show long correlation times of the order of forty seconds or 
Meeater while all other parameters have decay times of the 
Ku ir of two seconds. The single exception is the graph for 
5: “istor tvo during run five. The lead for this thermistor 
parted during the run and was replaced by the spare lead for 
5 last half of the run. This latter portion of the record 
was analyzed since it was the longer of the two records 
available. This proved to be a poor choice for there was 
intermittent contact to the circuit and therefore the record 
27 ce erratic. The consistency of the remaining records 
leads to the conclusion that the water was well mixed with 
respect to temperature for all runs. 

Generally an autocorrelogram shows fluctuations with the 
same kinds of periods as the original time series. The 
fluctuations have been put in phase so that they reach a 
Beximum at zero lag. Using the criterion that the spectrum 
of the time series probably has a general maximum near the 
first "peak" observed on the curve the following times are 
moted for the various parameters and runs as the general 
maximums. (For a complete listing of all peaks noted see 


fable 11). 


40 





I Run Three 


This run shows close agreement among the sound 
amplitude and thermistor one and thermistor two signals 
actually analyzed. The approximate time of the first peak 
of eleven seconds is valid for A and Tj whereas T, is 
slightly lower at about ten seconds. 

7 Lim Four 

There does not appear to be as close an agreement 
here as in run three. The first peak observed for A is at 
twenty-six seconds which is the same as the second period 
802”: the first at nine seconds. The salinity peak 
occurs at about 13.5 seconds while sound velocity exhibits a 
peak at 11.75 seconds. 

3. Run Five 

mis ls the most complete run available tor analysis 
and shows the most coherence among Signals. The first peak 
ranges from 7.61 seconds for the wave height to 7.08 seconds 
for the w component of turbulence. This close agreement 
27 “s is carried out through all five peaks plotted on the 
Graphs. 

nern Six 

Although this run contains a different set of 

parameters than run three it too exhibits the scatter shown 


mo run four. 


The conclusion drawn from this view of the analyzed data 
is that all parameters measured have the same general maxi- 


mam periodicity in the near surface environment. 
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C. POWER SPECTRUM COMPARISONS 

The use of the power spectra to show where most of the 
energy of the variance of the various parameters is produced 
agrees closely with the information derived from the auto- 


relation graphs. The information is found in Table III. 
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RUN THREE 

A HE 5 
GE 

T4 "E 

Тә  10* 

S 12.49 

RUN FOUR 

A None 

T> 9 

S None 

e ie. 7 5 

RUN FIVE 


AUTO-CORRELATION PEAKS 


(All times in seconds) 


LEDO 


20 


195 


26.24 


26 


13.54 


22.40 


Ei None evident 


A 7550 

T> 7.19 
* 

S 7 

u 7229 

v/ 7.08 

hp Dol 

с m 

BUN SIX 

m 83.54 

с ур 225 

5 16.08 


*Indicates peak not definite. 
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0.042 
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0.020 
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0.073 
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TABLE ILI 


POWER SPECTRUM PEAKS 


HZ 


HZ 


HZ 


HZ 


HZ 


HZ 


no peak 


02125 
0125 
0425 
0. 125 
0.125 


0.125 


0.062 
0.047 


0.088 


HZ 


HZ 


HZ 


HZ 


HZ 


HZ 


HZ 


HZ 


HZ 


23.8 Sec. Period 


50 Sec. Period 


16.1 Sec Period 
13.7 Sec Period 


12.0 Sec Period 


6.9 Sec Period 


8 Sec Period 


21.3 Sec Period 


11.4 Sec Period 
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TABLE IV 


CORRELATION DECAY TIMES 


Parameter Run Time to e | Correlation (sec) 
A 3 225 
Tj 3 >40 
Tə 3 z40 
S 3 Jus 
A 4 d 
To 4 740 
S 4 26 
C 4 1.9 
A 5 m2 
T > 38.4 
T> 5 lio. 36 
S 5 1.4 
u 5 176 
w 5 52 
ha 5 im 2 
C 5 14.7 
ET 6 740 
C 6 22 
S 6 419 
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Navy Undersea Research and Development Center 
Oceanographic Research Tower 


Figure l. 
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Carriage System Used 
to Mount Instruments 


Research Tower Schematic 


Figure 2. 
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Instrument Array 


Figure 8. 
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Thermistor and Velocimeter Mountings 


Figure 9. 
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